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Abstract. The d 2 d' configuration is analysed in group-theoretical terms. Starting 
from the table given by Condon and Odabasi for the configuration d 2 d! , we determine 
a set of convenient group-theoretical basis states, and rewrite the Coulomb matrix 
elements in terms of this new basis. Linear combinations from the different parts of 
the Coulomb operators are formed such that they have simple group transformation 
properties in our scheme. The sequence of groups that we use is f/(20) D SOt(3) x 
17(10) D 50 T (3) x SO s (3) x U(S) D SO T {3) x 50 s (3) x 50(5) D SO T (S) x SO s (S) x 
SOl(3), where T denotes the isospin of Simonis et al , in which electrons with the same 
angular momentum I but different principle quantum numbers n are accommodated by 
introducing the eigenvalue Mt of To. Using the Wigner-Eckart theorem and selection 
rules on the higher symmetry groups, the tables of the Coulomb and spin-orbit matrix 
elements for the reconstituted operators (with simple group transformation properties) 
are much simplified in terms of these basis states. 



1. Introduction 

The configuration d 2 d' has been observed in the excited spectra of some rather common 
atoms and ions. For instance, from the data of Sugar and Corliss (1985), Sc i, Ti n, 
V in and Cr iv all have excited levels with the configuration 3d 2 4d observed. In the 
case of vanadium and chromium, almost all the 27 possible spectroscopic terms of the 
configuration 3d 2 4d were found and their energies measured; for vanadium, a few 3d 2 5d 
levels were seen as well. In addition to the d 2 d' systems, the present work also applies 
to commonly found systems of the type d 8 d'. For example, 3d 8 4d has been observed in 
Co i and Ni n (Sugar and Corliss, 1985). For nickel, all 27 possible spectroscopic terms 
of the configuration 3d 8 4d have been found, as well as some levels of 3d 8 5d, 3d 8 6d, 3d 8 7d 
and 3d 8 8d. 

Condon and Odabasi (1980) calculated the Coulomb matrix elements of d 2 d' 
using classical methods. Their basis couples the two equivalent electrons (d 2 ) to an 
intermediate state ^L', to which the d' electron is coupled to give the final state ^L. 
The basis is denoted by \(^L')^L), where [S] = 25" + 1, the multiplicity. This basis is 
quite good, as the d' electron usually interacts rather weakly with the other d electrons. 
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As one can check from the data for V in and Cr iv (Sugar and Corliss, 1985), there is 
not a substantial mixing of the states in this basis except for 2 P, 2 D and 2 F. 

Instead of this physical basis, we choose to use a different one that has simple 
transformation properties under higher symmetry groups. There are advantages in doing 
so, particularly when they are taken as intermediate states in a perturbation expansion 
(see, for example, Hansen et al 1997). In that case, the sum over intermediate states that 
appear render the choice of basis irrelevant to the physical properties being investigated. 
And certainly, if we introduce higher symmetry groups into the analysis, even though 
they might not give "good" quantum numbers, in the sense that they do not commute 
with the full Hamiltonian, they give simple results as far as the computation of matrix 
elements is concerned. The selection rules, and especially the Wigner-Eckart theorem, 
when applied to higher groups, strikingly reduce the amount of calculation and expose 
new mathematical structure. 

Other than matrix elements among the states of the d 2 d' configuration itself, our 
interest in configuration interaction (CI) means that we will also need the matrix 
elements that connect d 2 d! to d 3 . The physical basis that is commonly used for d 3 
coincides with the group-theoretical basis that utilizes the higher unitary and symplectic 
groups C/ (10) and Sp(10). For this reason, it is best to use the same group-theoretical 
basis for the configuration d 2 d' as that needed to study the configuration interaction 
between d 3 and d 2 d'. Once again, the Wigner-Eckart theorem and selection rules for the 
higher symmetry groups are invaluable. 

The isospin formalism (Simonis et al 1984, Kaniauskas et al 1987) that we will 
introduce into our work also provides a very useful setting for the study of CI. In the 
isospin formalism, the states in d 3 and d 2 d' that carry the same group labels belong 
to the same isospin multiplet, but with different isospin z-component M T . We can 
therefore easily obtain the CI matrix elements from the ones for d 2 d' if everything, the 
states and the operators, have simple transformation properties under the symmetry 
groups. These considerations give us strong reasons to abandon the physical basis for 
d 2 d' and turn to the group-theoretical one instead. 

2. The basis states 

The classification used in Condon and Odabasi (1980), in group-theoretical language, is 

SO d s {3) x SO d L (3) x SO d s(3) x SOf (3) D SO s (3) x SO L {3). (1) 

As explained earlier, this provides a set of good physical states, but there are also good 
reasons to use an alternative basis which carries a set of simple group labels, that is, 
irreducible representations (irreps). The groups that we have chosen are 

17(10) D SO s {3) x U(5) D SO s {3) x SO(5) D SO s {3) x SO L {3). (2) 

The groups U (5) and 5*0(5) are based on the five orbital states of a d (or d') electron; the 
group £7(10) utilizes the spin of the electron as well. Another more common alternative 
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New basis 



Old basis in terms of \d 2 ( 2S ' +1 L'),d'; 2S+1 L) 



| 2 [21][3](30) 2 S 
| 4 [111][21](21) 2 P 
| 2 [21][111](11) 2 P 
| 2 [21][21](21) 2 P 
| 4 [111][111](11) 4 P 
| 2 [21][21](21) 4 P 
| 4 [111][21](10) 2 P 
| 4 [111][21](21) 2 P 
| 2 [21][21](21) 2 D 
| 2 [21][21](10) 2 P 
| 2 [21][3](10) 2 P 
| 2 [21][21](10) 4 P 
| 2 [21][21](21) 4 D 
| 4 [111][21](21) 2 P 
| 2 [21][111](11) 2 P 
| 2 [21][21](21) 2 P 
| 2 [21][3](30) 2 P 
| 4 [111][111](11) 4 P 
| 2 [21][21](21) 4 P 
| 4 [111][21](21) 2 G 
| 2 [21][21](21) 2 G 
| 2 [21][3](30) 2 G 
| 2 [21][21](21) 4 G 
| 4 [111][21](21) 2 P 
| 2 [21][21](21) 2 P 
| 2 [21][21](21) 4 P 
| 2 [21][3](30) 2 / 

| 2 [21](21)(10) 2 D 
| 2 [21](1}(10) 2 £> 
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For states belonging to [111] 10, we use the same phases as in d 3 ; otherwise, we make 
an arbitrary phase choice when there is no precedent to guide us. 



used in configurations of equivalent electrons is the sequence 

U(10) D Sp(10) D SO s {3) x SO{5) D SO s (3) x SO L {3) (3) 

where the symplectic group Sp(iO) is closely related to the seniority scheme widely used 
in atomic spectroscopy. It is, however, felt that the idea of seniority might not as useful 
in configurations with inequivalent electrons. At least not in the way we used it in 
configurations of equivalent electrons. 

To simplify the notation, we use square brackets for unitary group labels, with a 
subscript to denote the order of the group. The zeros in the labels will either be omitted 
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or abbreviated as a dash (-). Also, we write the negative values with an overline to 
save space. For example, the irreps [210 s ], [1 2 6 - l 2 ] of U(10) and [11100] of £7(5) will 
be denoted as [21]io, [11 — TT]io an d [Hljs respectively. The exponent form is retained 
if the same number (other than zero) occurs more than three times. The advantages 
of using such a notation will become clear when we discuss the generic branching rules 
of the type U(2n) -> SU{2) x U(n) = SO (3) x U(n) in table |; as well as to tabulate 
Kronecker products for unitary groups of any order n (see table 

The mapping of states from the classical basis to the new basis was worked out by 
Judd (1997) and the author, and is given in table |l|. The reader should for the moment 
disregard the first superscript labeling each state. That is the isospin label which will be 
introduced later in section f|. If the alternative scheme involving the symplectic group 
Sp(10) is used instead, almost all of the states in the first column of table |l| corresponds 
to a single state in that scheme. The only mixing comes from a pair of 2 D states that 
share the same labels in U(10) and SOs(3) x SO(5). They are given on the last two 
lines in table [|. 

3. The Coulomb interaction matrix elements 

The matrix element table given in Condon and Odabasi (1980)|§] can be transformed 
to the new basis without difficulty At the same time, we also followed Judd (1998) 
and Racah (1954) and take specific linear combinations of the two-electron Coulomb 
operators so that they transform irreducibly under the action of the group 50(5). 

The full two-body Coulomb interaction contains three distinct parts. One involves 
only the <i-electrons (or only the cf-electrons, if we had more than one cf -electron.) The 
other two are the direct and exchange parts that involve one d and one d! electron. 

We can write the perturbing Hamiltonian as 

Hi = J2 /*( d » d)F k (d, d) + f k (d, d')F k (d, d!) + g k (d, d')G k (d, d') 

k 

= £ + J.T + gp l (4) 

i 

where the sums run over k = 0,2,4 and % = 0,1,2. The F k and G k are the usual 
direct and exchange radial integrals, and the matrix elements of fk{d, d') and g k (d, d') in 
d 2 d' in the classical basis are tabulated by Condon and Odabasi (1980). The operators 
e , e x and e 2 on the second line, constructed from the operators f k (d, d) by Judd (1998) 
and Racah (1954) independently, transform irreducibly as (00), (00) and (22) of 5*0(5) 
respectively. They are given by 

e = f (d,d) 

e l = 7 -[f (d,d) + 5f 2 (d,d) + 9U(d,d)] 

e 2 = 63[f 2 {d,d)-U(d,d)\ (5) 

§ A negative sign for the matrix element (( 3 : F) 2 H\go\(} G) 2 H) is missing in their table; its correct value 
should be -105. 
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and the corresponding E 1 are 

7 fiS 
E° = F (d, d) - -F 2 (d, d) - -F 4 (d, d) 

5 4^ 
E 1 = -F 2 (d,d) + —F A (d,d) 

E 2 = ~F 2 (d,d)-~F 4 (d,d) (6) 

where F = F°, F 2 = F 2 /49 and F A = F 4 /441 as defined by Condon and Shortley (1953). 
Replacing fk(d, d) by fk(d, d') and gk(d, df) in equation (^), we arrive at the operators f i 
and ^ respectively; and they transform as (00) and (22) in SO (5) just like the operators 
ej. Similar substitution in equation @ gives us the relation between F and the direct 
radial integrals F k (d,d'); and likewise, the relation between (? and G k (d,d'). 

In terms of the creation and annihilation operators, or rather in terms of 
v^J = \/2w^ = —\^2(d)dY oh \ the operators ej are given by 

eo = 2 : w £t ■ v £ } : ( 7ffl ) 
Pl _ I . „(o) . (o) . , . (2) . (2) . . (4) (4) . 

e 2 = 9:^-^: -5 : v$ ■ v$ : (7c) 

where colons denote normal ordering!]]. For instance, with Einstein summation 
convention over repeated Greek indices (/^, v etc.) used throughout the paper, the 
operator eo from the above can be written as 

eo = \ ■ vf] ■ vf] : = - : (dJ,d M ) • (d\d v ) : = — {S^d^d,,) = ^(44^) • ( 8a ) 

The direct and exchange operators, j i and g i} also admit similar expressions in 
terms of the operators v^(= —^(d^d'Y ^), and v$ d ,. For example, / and g Q 
are given as 

fo = 5:v®.v$l: = (d^d'M (8b) 

3o=5:«S-«S°J: =(d\d^dj p ). (8c) 

The extra factor ~ in equations (|7aD and ( |8c| ) is to allow for the pairwise interaction 
between the identical <i-electrons being counted twice. The other direct and exchange 
operators f i and 'g i are similarly given as in equations (|7^) and flTcl), with an extra factor 
of two. As a note, normal ordering removes the self- interaction terms in ef, and it take 
cares of the ambiguity that arises in the case of the exchange operators 'g i , in which 
v dd> " v d/d an d v d>d ' v dd> are both legitimate but distinct forms, if not normal ordered. 
In terms of the new basis states given in table [IJ, the matrix elements of the operators 
f i and g { can be easily obtained from the table by Condon and Odabasi (1980). We 

|| The colon notation is widely used in Quantum Field theory; see Weinberg (1995) for example. 
Lindgren and Morrison (1981) also used normal ordering in atomic theory, and denote it by curly 
brackets. We prefer colons to avoid confusion with anti-commutation brackets and brackets in general. 
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(11) 4 F 


1 


2 


2 


2 [21] [21] 


(21) 4 F 


1 


2 


-1 


4 [111][21] 


'10) 2 L> 


1 


2 


2 


4 [111][21] 


'21) 2 L> 


1 


2 


2 


2 [21][21] 


'21) 2 L> 


1 


2 


-1 


2 [21][21] 


'10) 2 L> 


1 


2 


-1 


2 [21] [3] 


10) 2 L> 


1 


2 


-1 


2 [21] [21] 


'10) 4 L> 


1 


2 


-1 


2 [21] [21] 


'21) 4 L> 


1 


2 


-1 


4 [111][21] 


(21) 2 F 


1 


2 


2 


2 [21][111] 


(11) 2 F 


1 


2 


-1 


2 [21] [21] 


(21) 2 F 


1 


2 


-1 


2 [21][3] 


(30) 2 F 


1 


2 


-1 


4 [111][111] 


(11) 4 F 


1 


2 


2 


2 [21] [21] 


(21) 4 F 


1 


2 


-1 


4 [111][21] 


(21) 2 G 


1 


2 


2 


2 [21] [21] 


(21) 2 G 


1 


2 


-1 


2 [21][3] 


(30) 2 G 


1 


2 


-1 


2 [21] [21] 


:2i) 4 G 


1 


2 


-1 


4 [1H][21]( 


2l) 2 H 


1 


2 


2 


2 [21][21]( 


21) 2 H 


1 


2 


-1 


2 [21][21]( 


21) 4 H 


1 


2 


-1 


2 [21][3] 


(30) 2 / 


1 


2 


-1 



7 
3 



7 

3 
VI 
3 



2 

1 1 



1 1 





1 1 
1 1 






7 

3 




7 

3 
VI 
3 



10 10 

10 10 
2 


1 1 
1 1 
2 


1 1 
1 1 


2 



VI 

3 




VI 

3 
13 

3 



14 

3 




_ 7 
3 

V7 
3 



4 

2 0-1 

-10 2 



3 





-1 
2 


- 






14 
3 

VI 

3 



7 
3 
2 0-10 




-10 2 


3 
2-10 
-12 
4 
3 

2 -1 
-1 2 
3 
4 



VI 

3 




VI 

3 
26 
3 



-2 
2 0-1 

-1 -1 

-3 



14 
3 




_ 7 

3 

2^7 
3 



o -1 



o 4 

-7 
-3 
2 0-10 

-10-1 
-2 

-3 
2-10 
-1 -1 
0-2 
-3 
2 -1 
-1 -1 
-3 
-2 



_ VI' 

3 





VI 

3 

_ 13 
3 . 



The matrix elements go from Condon and Odabasi's table possess a variety of values. 
In the new basis, the corresponding ~g Q matrix elements are diagonal and only take the 
values 2 and —1. 



simply transform their matrix elements^ of fk{d, d), fk{d, d') and gk{d, d') into the new 
basis, then take the corresponding linear combinations as given in equation In this 
way, we arrive at table |2j. 

3.1. The irreps of the Coulomb operators 

As can be seen from table [l], all the states have definite group labels with respect to 
the symmetry groups in (0). We have now to determine the transformation properties 
of the nine operators e,, f i and Of course, we know that they are scalars in the 
spin and orbital spaces. As mentioned in the previous section, these nine operators thus 

% Notice that the matrix elements of f k (d,d) in the table of Condon and Odabasi (1980) are simply 
taken from the configuration d 2 . 
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" -2 4/7 -7 
4/7 7 -4/7 
. -7 -4/7 -2 
7 -4/14 
-4/14 5 
2/21 -V21 
2/21 4 -8 -/21 -2/3 
-/21 -8 4 2/21 -2/3 
-/21 2/21 
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-7 -6/2 8 2/5 
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" -16 11 2/55' 
11 -16 2/55 
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-9 
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. 7 - 2 
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5 



-18 
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16 -6/2 7 -2/5 
-6/2 -6 6/2 
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-11 -32 -2/55 
-2/55 -2/55 -14 
-7 
4 -7 
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-4 -4/7 7 
4/7 -7 -8/7 
7 -8/7 2 
14 4/14" 
4/14 16 
4/21 /21 " 
4/21 8 8 /21 2/3 
/21 8 -4 -2/21 -4/3 
/21 -2/21 
2/3 -4/3 
" 
12 

16 -6/2 7 -2/5" 
-6/2 3 -12/2 

7 -12/2 -8 4/5 
-2/5 4/5 5 
-6 12" 
12 6 

-32 -11 -2/55" 
-11 16 4^/55 
-2/55 4/55 7 
-2 
-4 -71 
-7 2 
-12 
-5 



constructed belong to either (00) or (22) of 5*0(5). We will now find the {7(5) and £7(10) 
labels of these operators. 

The two-body operators contain two creation and two annihilation operators. As 
the creation operators {dft or w) transform as while annihilation operators [T]s » 
the product 

[l]e x [1] 5 x [T] 5 x [T] 5 = 2[0] 5 + 4[1 - Y] 5 + [2 - TT] 5 + [11 - 2] 5 + [11 - TT] 5 + [2 - 2] 5 

must contain the appropriate labels for the two-body operators. We only need to pick 
out those t/(5) irreps from the Kronecker product above that contain either (00) or (22) 
of 5*0(5). The branching rules {7(5) — > 50(5) are given by Wybourne (1970)Q, for 
instance. We are then left with [0]s, [11 — ll]s and [2 — 2] 5 where the first irrep contain 

+ The table by Wybourne (1970) does not include U(5) irreps with negative values. But if we look at 
the reduction U(b) — > SU(5) — > SO(5) instead, the fact that [ax, a%, ■ ■ •] and \a\ + c, + c, • • •] of U(5) 
reduce to the same SU(5) multiplet tells us they must contain the same 5*0(5) components. So, we 
can easily remove the negative numbers in each irrep of U(5) without changing its SO(5) content. 
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U{2n) 50(3) x U(n) 



[Ihn 2 [l]n 

[H] 2 „ 3 {U] n + 1 [2] n 

[l r hn '•+ 1 [11, i + r - 1 [21 ( ^ 2) ]« + --- 

[2] 2 „ 1 [H]n+ 3 [2]„ 

[21] 2n 2 [lll]„ + 2 ' 4 [21] n + 2 [3] n 

[1-T] 2 „ 3 [0]„ + 1 - 3 [l-T]„ _ 

[11 - TTj 2n 1 - 5 [0]„ + ; • 1 - l]„ + i> 3 > 5 [11 - TT]„ + 3 [ll - 2]„ + 3 [2 - TT]„ + - 2] 

[2 - 2] 2n LB [0]„ + L ' 3 ' 3 ' 5 [1 - T]„ + 1 [1 1 - TT]n + 3 [11 - 2]„ + 3 [2 - TT]„ + !- 3 ' 5 [2 - 2] 

[2 - TT] 2n 3 [0]„ + ^^[l - T]„ + 3 [11 - TT]„ + 1 [11 - 2]„ + ^ 3 - 5 [2 - TT]„ + 3 [2 - 2\ n 



(00) of 50(5), the second one contain (22) and the last one both. 

On the U(10) level, we use the branching rules given in table |3] and pick out the 
irreps in the Kronecker product [1] 10 x [1] 10 x [1] 10 x [T] 10 that contain either 1 [0]5, 

— TT]5 or x [2 — 2] 5 of SOs(3) x £7(5). (The Coulomb operators are spin scalars.) 
The eligible ones are [0]io, [11 — TT]io an d [2 — 2] 10 . 

The group labels for the nine operators e i: f { and 'g i can thus be written down 
rapidly as given in table f|. (Once again, the reader can ignore the isospin superscript 
for the time being.) The only ones that need explanation are the three e« operators as 
well as f and ~g Q . The three equivalent electron operators e« do not contain the irrep 
[2 — 2] io of £7(10) because the two creation operators (and likewise, the annihilation 
operators) must form an antisymmetric product in £7(10). So no [2 — 2] 10 appears. In 
fact, they were examined by Judd and Leavitt (1986), and it is found that eo belongs 
to [0]io of £7(10), e 2 belongs to [11 — H] 10 , and e\ is a mixture of [0]i and [11 — TT] 10 . 

For the operators eo, f and g Q , they can be shown to be £7(10) scalars rather 
easily. From the expressions (|8^) through (|3c]), it is a simple exercise to show that they 
commute with all the 100 generators d„ + d'^d' u of the group £7(10). (See section £| for 
further discussion of the group generators.) Hence, they are all £7(10) scalars. 

3.2. Selection rules and the Wigner-E chart Theorem 

With the group labels for the operators worked out, we are now in a position to 
appreciate the use of the new group-theoretical basis. We can look at the situation 
on different group levels. First of all, the six 5*0(5) scalars operators eo,ei, f , f x , 
~g and 7^ cannot connect states with different 5*0(5) labels. And for the other three 
operators e 2 , f 2 and g 2 which belong to (22), the Kronecker products for 50(5) account 
for all the vanishing matrix elements one can find in table |2| on those three columns. 

When we apply the Wigner-Eckart theorem on the group £7(10), we expect a lot of 
simple proportionality relations among the different operators. Between the operators 
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Table 4. Group labels of the Coulomb interaction operators. 



operator 


50t(3) x E/(10) 


t/(5) 








50(5) 


eo, e 


x ' 8 [0]io 


[0] 5 








(00) 


ex, ei 




[0] 5 








(00) 




^[ll-ITjio 


[0] 5 + [2 - 


^2] 5 




(00) 


e 2 , e 2 


^ 5 [ll-TT] 10 


[11 - 


TTk 

JO 


+ [2 


~2] 5 


(22) 


foi 9o 












(00) 


fx: 9x 


^[Ojio 


rou 

L J ° 








(00) 




^[ll-TTlxo 


5 H 


- [2- 


-2] 5 




(00) 




1 [2-"5] 1 o 


5 H 


- [2- 

L 


-2] 5 




(00) 


fit 92 


^[H-TTlio 


[11- 


TT] 5 


+ [2 


~2] 5 


(22) 




M2-2]io 


[11- 


H] 5 


+ [2 


-2] 5 


(22) 


e 


Moiio 


[0] 5 








(00) 


ex 


Moiio 


[0] 5 








(00) 




x [2-2]io 


[0]s H 


-[2- 


-2] 5 




(00) 




1 f9 91 

— zj io 


[11- 


11] 5 


+ [2 


91 


(22) 


eo + eo 


Moiio 


[0]5 








(00) 


ex + ex 


Moiio 


[0] 5 








(00) 




Mn-T^io 


[0] 5 H 


-[2- 


-2] s 




(00) 


e 2 + e 2 


Mll-T^io 


[11- 


TT] 5 


+ [2 


-2] 5 


(22) 


2e - e 


5 [0]io 


[0] 5 








(00) 


2ei — ei 


5 [0]io 


[0] 5 








(00) 




5 [ll-TTjio 


[0] 5 H 


-[2- 


-2] s 




(00) 


2e 2 - e 2 


5 [H-TT]io 


[11- 


TT] 5 


+ [2 


-2] 5 


(22) 



f i and g i7 we found 

([lll]io/3|7il[Hl]io7) = ([lll]io/3|^|[Hl]io7) (9a) 

<[lll] 10 y9|7<|[21]ioT) = <[lll]io/?|&|[21]io7> (96) 

for i = 0,1 and 2; and (3, 7 denotes the rest of the quantum numbers needed to specify 
the state. As for the relations with the operators ej, we have 

([lll] 10 /3|7il[lll]io7> =2 ([lll]io/3|e i |[lll] 10 7) (10a) 

([lll] 10 /3|7,I[21] a oT> = -<[lll]ao/3|e<|[21] a0 T> (106) 

([21] 10 /3|7i + ^|[21]io7> = ([21] 10 /3| ei |[21] 10 7) • (10c) 



Most of the proportionality relations above are the result of the Wigner-Eckart 
theorem at work on the group {7(10). To demonstrate, let us use a few examples. 
For % = 2 in equation (|9a|), since [111] 10 x [HI] 10 contains [11 — 11] 10 once and no 
[2 — 2] 10 appears (see table |]), only the [11 — IT] 10 part (but not the [2 — 2] 10 part) of the 
operators f 2 and g 2 contributes to their matrix elements. Hence, by the Wigner-Eckart 
theorem, they should be proportional to each other, and the proportionality constant 
turns out to be unity. Analogous argument applies to explain the proportionality 
relations in (|91|), ( |10a| ) and ( |106| ) for i = 0,2. The i = 1 operators contain the irreps 
[0]io and [11 — IT] 10, which both contribute to the matrix elements in general. So, the 
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Wigner-Eckart theorem fails to provide an explanation for these i — 1 cases in ( |9^ ) 
and ( |10g| ); although it works well with (^) and ( |10^j ) where the scalar part [0]i does 
not contribute to the matrix elements. But in any case, it cannot explain why the 
proportionality constants are the same for i = 0, 1 and 2. 

The last relation ( |10c] ) is a little bit different. Since [2T]io x [21] io contain both 
[11 — ll]io and [2 — 2] 10 , we do not have a simple relation as we had before. For i = 2, 
if we write equation (|10 c|) as (g 2 ) = (e2) — (/ 2 )> we can interpret (62) and (f 2 ) as two 
independent sets of isoscalar factors; hence (g 2 ) can be written as a linear combinations 
of the two sets of isoscalars. Once again, the argument does not apply to the i = 1 case; 
and for i = 0, the three sets of matrix elements (eo), (f$) and (g ) are proportional to 
each other as the operators are all U(10) scalars. Despite these differences between the 
i = 0, 1 and 2 cases, they display the same relation as given in (|10c] ). We will explain the 
above relations in a more elegant way in the next few sections. Relations (0) and (|9|) 
will be explained in section |373"| , and relations (|10a|) - (|10dj in section [O] making use of 
the isospin structure. 

Before we go on any further, let us pause and reflect on what we have achieved 
so far. We are trying to find all the matrix elements of the nine operators e«, f { and 
Tji for the configuration d 2 d' . The matrix elements (e^) are rather well known, as they 
come strictly from the d 2 matrix elements (see Condon and Odabasi, 1980). Now with 



equations (0), (0) and ( |10dJ , the g i matrix elements can be easily obtained from 
those of f i and e^. Furthermore, with equations ( |10a| ) and (|10£| ), the matrix elements 
([111] lo/S I fi\ [HI] 107) and ([lll]io/3|/j| [21] 107) are readily obtained from the matrix 
elements; we are only left with ([21] 10 /?|/j| [21] 10 7) to work on. Is there an easy way 
to find these matrix elements? By observation, we find that the matrix elements of 
e-i = Ci + fi are surprisingly simple. For % = 0, the sum e + f is always 3; for % = 1, 
the sum e\ + f 1 is diagonal, and the value depends only on the U(5) and 50(5) irreps, 
but not the £7(10) label nor the spin and orbital ranks. For the case i = 2, the sum is 
almost diagonal, with the three exceptions 

(d 2 d' [T] U[21]{21) [s] D\e 2 + J 2 \d 2 d' [T] U[21]{10) [S] D) = 6^21 . 

The diagonal values once again do not depend on the U(10) irrep and the spin, although 
they do depend on the orbital rank. These properties, which will be addressed in 
section greatly simplify our task of finding the f i matrix elements. With these 
simple relations, (almost) all the matrix elements can be related to the known matrix 
elements. We will elaborate this further in section || 



3.3. The operators and with simpler U(10) transformation properties 

As mentioned before, the ten creation operators d^ belong to the irrep [l]i of £7(10); so 
do the ten components d'^. When we form products from the two independent sets, we 
can get [2]io or [ll]io- From the knowledge of unitary groups and their representations, 
we know that [2] 10 corresponds to the symmetric product, while [11] 10 the antisymmetric 
one. More precisely, (c^d'J + d'^d\) and [d^d^ + d'^d u ) have 45 components each; they 
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are the antisymmetric products and belong [11] 10 and [TT] 10 respectively. The other 
combination (d^ d'J — d'^d) u ) belongs to [2]i , has 55 components; and (d^d' u — d' d u ) 
belongs to [2]io- Of course, the symmetric product [2] 10 is identically zero for two 
equivalent electrons. 

In terms of creation and annihilation operators, /, and g { can be written as 

Ji = Z 2a ^ ■ v dd ■ v d% ■ ~ C^C^d^d^dy 

k 

9i=Yl 2a ^ ■ v ddi ■ v d/d ■ ~ Z! C^C^d'^d^ (11) 

k 

where are the corresponding coefficients as given in equations (|7o|) to 
are the appropriate Clebsch-Gordan (CG) coefficients. They are invariant under the 
simultaneous interchange of d^ with dT and d with d'. But they are neither symmetric 
nor antisymmetric with respect to the interchange of the creation operators alone. In 
fact, if we interchange d^ and rf'f (or d and d') alone, f { becomes g i: and vice versa. If 
we take the combination 

e» = ^(7. + 9%) ~ J2 C^C^dld'^(d' ( d u + d c d' u ) 

we first note that the annihilation operators within the parenthesis transform as the 
antisymmetric product [TT]io. With a little algebraic manipulation, we can obtain a 
similar expression with the creation operators grouped as a [ll]io product. This tells 
us the new operators ej belong to [11 — ll]io of {7(10). Similarly, ej = \{fi — ~g~i) is the 
symmetric product, and belongs to [2 — 2] 10 . There are flaws in the above argument; 
namely, we ignored the scalar parts. The arguments above do not exclude the possibility 
that the operators can be {7(10) scalars. As we already know, f and g are both {7(10) 
scalars, so the sum or difference (e or e ) must also be {7(10) scalars; although their 
expressions look as if they belong to [11 — TT]io and [2 — 2]i respectively. Similarly, e\ 
and e"i can contain a scalar part and, in fact, they do. The transformation properties 
of these new operators and e, are summarized in table £|. In terms of these six new 
operators, the perturbative Hamiltonian in (|J) becomes 

H x = Y^E' + J.T + g^ 

i 

= J2 e i Ei + ^^(T + (?) + li^k(T - cf) 

= Y f e i E i + e i E i + e i E i . (12) 

i 

From a group-theoretical point of view, it is more convenient to use these new 
operators, and e^, with simpler transformation properties, rather than the distinct 
direct and exchange operators, f i and g i . Notice that the operators have the same 
transformation properties as ej. The matrix elements of e-i and e"j are given in table 
in which we copied down the ej matrix elements from table [2] to make it complete. 
These new operators transform irreducibly in {7(10), makes applying the Wigner-Eckart 
theorem as well as the selection rules on the {7(10) level so much easier and more elegant. 



7c), and C 



in' 
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Table 5. Matrix elements for the operators e,-, e^, and 



e e e e 



ei 



ei 



ei 



ei 



2 [21][3](30) 2 S 
4 [111][21](21) 2 F 
2 [21][111](11) 2 F 
2 [21][21](21) 2 F 
4 [111][111](11) 4 F 
2 [21][21](21) 4 F 
4 [111][21](10) 2 L» 
4 [111][21](21) 2 J D 
2 [21][21](21) 2 L> 
2 [21][21](10) 2 L> 
2 [21][3](10) 2 J D 
2 [21][21](10) 4 L> 
2 [21][21](21) 4 L> 
4 [111][21](21) 2 J F 
2 [21][111](11) 2 J F 
2 [21] [21] (21)^ 
2 [21] [3] (30)^ 
4 [lll][lll](ll)fF 
2 [21] [21] (21) 4 F 
4 [111][21](21) 2 G 
2 [21][21](21) 2 G 
2 [21][3](30) 2 G 
2 [21][21](21) 4 G 
4 [lll][21](21) 2 i? 
2 [21][21](21) 2 # 
2 [21][21](21) 4 H 
2 [21][3](30) 2 / 



12 3 



3 
3 
3 
3 
3 
3 
3 
3 
3 
3 

3 o 
2 

3 



7 

3 




7 

3 

2^/7 
3 



2 

1 1 



1 1 





1 1 
1 1 






7 

3 




7 
3 
V7 

3 



10 10 

10 10 
2 


1 1 
1 1 
2 


1 1 
1 1 



2 



2y/7 

3 




2y/7 

3 

13 
3 



14 

;s 




_ 7 
3 

2y/7 

3 



1 

2 0-1 

-10 i 






-1 



1 



7 

6 

VI 

3 



1 



1 







2 0- 

-10 \ 



2-10 
-1 I 

o 6 i 



2 -1 

-1 



1 



2 J 



2\/7 
3 




£L 

3 

13 
6 






oo I 





3 




oof 


3 




7 

2 



,0 -/7 
7 
3 
0" 

| 
3 

3 

| 
3 
3 


^ 

u 2 

3 
3 



" 




13 

2 - 



6 
3 

3 

3 
7 
3 
3 
7 

13 
7 
3 
3 

3 
G 

3 
3 
3 
G 
3 
3 
3 
3 
G 



For example, equations ( |9q ) and ([JtJ) are now written as 

([lll]io/3|ei|[lll]io7> = (13a) 
([lll]io>9|ei|[21] loT ) =0. (136) 



From the Kronecker products given in table |6|, neither [llljxo x [Hl]io nor [Hljio x [21] io 
contains [2 — 2] 10 . Therefore, equation ( 13£|) as well as the % = 2 case in ( p.3o| ) come from 
the U(10) group selection rule. The matrix elements (So) = in (|13a|) can be obtained 



using the derived eigenvalue expression for eo in section |4.6| after we introduce the isospin 
structure. And (ei) = in ( |13q| ) can be verified using the results in section [5]. Most of 
the Wigner-Eckart relations we had before are now re-interpreted as the selection rule 
on the group {7(10). This also explains why the proportionality constants are the same 
for all i, and why they are 1. 

The other three new operators have identical group labels with the corresponding 



Coulomb and spin- orbit interaction matrix elements in d 2 a" configuration 
Table 5. laterally continued. 



13 



C2 



e-2 



(■2 



-9 

" -2 V7 -7 
V7 7 -V7 
. -7 -V7 -2 
7 -V14 
-V14 5 
2^/21 -y/21 
2^/21 4 -8 -y/21 -2^3 
— v/21 -8 4 2-J/21 -2^3 
-721 2^21 

o -273-273 

3^21" 
3^/21 12 
8 6^2 -7 275" 
6^2 -3 -6^/2 
-7 -6y/2 8 2^/5 
2^5 2^5 -5 . 
-3 -12" 
-12 15 
-16 11 2^/55" 
11 -16 2^55 
2V55 2V55 -7 . 
-9 

-2 7 " 
7 -2 
-9 
5 



_9 
2 

-4 -477 7 
-477 \ -2^7 
7 -277 -1 
14 4714^ 
4714 f 

4721 721 

4721 8 8 721 273 
721 8 2 721 -73 

721 721 
273 -73 
n V 21 " 

u 2 
V21 „ 

2 u 

16 -672 7 
-672 -§ -372 
7 -372 4 

-275 75 

-6 12 



-275- 




12 f 



-32 



-11 -2755 



-11 -8 
-2755 755 

_9 
2 

-4 -7 
-7 -1 

_9 
2 

5 
2 



755 

_ 7 
2 



_ 27 
2 



f 677 

677 -3 




-21 
u 2 . 

1 
I 







3721 
00 373 



6 3721 373 







2 





3y/21 

2 

-6 










-375 
-375 -f 





-f 972 

972 12 



l 


-24 

-3755 

_5 
2 



-3 

15 
2 

15 
2 





-3755 

_ 21 
2 



-27 

-6 

21 

-6 

21 

-6 

6721' 

6721 12 
12 6721" 

6721 



6721" 

6721 12 
24 

-9 

24 
-15 

-9 

24 
-16 
-16 

-7 
-16 

-6 

-6 

-6 
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ej operators. Using the Wigner-Eckart theorem, we find 

([lll] 10 /3|e<|[lll]ioT) = 2([lll] 10 /9|e i |[lll]io7> (14a) 

([lll]io/3|e<|[21]io7) = " ([lll]io/3N[21] 10 7> ( 14 ^) 

(filUPlemUl) = i([21] 10 /3|e,|[21] lo7 ) (14c) 



which are the same as equations ( lOcj) to (|10d) . Once again, the i = 1 cases in Q14q|) 



and (|14tj) do not come from the Wigner-Eckart theorem. In a sense, we are just rewriting 
those three equations without adding much understanding to the problem we are facing. 
But introducing these new operators and e"j is a very crucial step towards the full 
understanding of the underlying group structure to the problem on hand. 

4. Isospin and the group U(20) 

Simonis et al (1984) first introduced the idea of isospin into the theory of atomic 
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Table 6. Some Kronecker products for Unitary Groups U(n). 



[ii x rn 


[21 + [111 






fll x [T| 


[1 - T] + [0] 






[111 X [111 


[221 + 12111 + [1 4 1 

L J 1 L J 1 L J 






[in x im 

L J L J 


[11 -TT] + [1 -T] + [0] 






[111 x [21 


[311 + [2111 

L J 1 L J 






[111 x [21 

L J L J 


fii — 21 + [i - n 

L J 1 L J 






[21 x [21 

L J L J 


[41 + [311 + [221 

L J 1 L J 1 L J 






[21 x [21 


[2 - 21 + [1 - Tl + [01 

L J 1 L J 1 L^J 






run x run 

L- L - L - L J [■ i -- L - L J 


[2221 + 122111 + [21 4 1 + [1 6 1 

L J 1 [ « J i [ J 1 L J 






fllllxfTTTl 

L J L J 


[ill - TTT] + [11 - TT] + [1 - T] + [0] 






[111] x [21] 


[321] + [3111] + [2211] + [21 4 ] 






[111] x [21] 


[111 - T2] + [11 - 2] + [11 - TT] + [1 - 


1] 




[lll]x[3] 


[411] + [3111] 






[lll]x[3] 


[111-3] + [11 — 2] 






[21]x[21] 


[42] + [411] + [33] + 2[321] + [3111] + 


[222] + 


■ [2211] 


[21]x[2T] 


[21 - 15] + [11 - 2] + [2 - TT] + [2 - 2] 


+ [H- 


-TT] + 2[1-T] + [0] 


[21] x [3] 


[51] + [42] + [411] + [321] 






[21] x [3] 


[21 - 3] + [11 - 2] + [2 - 2] + [1 - T] 






[3]x[3] 


[6] + [51] + [42] + [33] 






[3]x[3] 


[3 - 3] + [2 - 2] + [1 - T] + [0] 







The above products are obtained using the method of Young's tableaux. When 
negative integers are involved in the irrep, we first add a constant to it so that all 
integers becomes non-negative; then use the standard Young's method to find the 
Kronecker product, and finally subtract the same constant from the resulting irreps. 
The Kronecker products are valid for large enough n. For instance, [111 — lll]g 
contains at least six entries and is identically zero for U(5). 

spectroscopy. The idea, borrowed from nuclear spectroscopy, is to treat two electrons 
with different principal quantum numbers n as two states (f or j) of a generic (^-electron 
in the isospin space. This gives an elegant explanation to the simple relations observed in 
the previous sections, as well as providing new results such as the closed-form expressions 
for the scalar operators. Furthermore, putting d and d' electrons on the same footing 
lets us conveniently switch between configurations d 3 , d 2 d' and so forth; hence CI can 
be taken into account in a natural manner. 

4-1. Motivation 

In an attempt to exploit the orthogonality relation (Judd and Leavitt, 1986 and Judd 
et al 1982) between the nine operators we are looking at (either the set e iy f { , g { or e iy 
Ci and ei), it becomes clear that the situation we are now facing is very different from 
the equivalent-electron case. In the latter one, the operators in question (e$) belong 
solely to a single irrep [11 — TT\ W of U(10), with the exception of a complete scalar 
operator, eo (see table £|). The three operators eo,e[ — e\ — |eo and e2 are orthogonal 
to each other. Orthogonality requires summing over all states with the irrep [l^jio (all 
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the states in d 3 , for instance). But now, the states in d 2 d' belong to two different irreps 
[llljio and [21] io- Another very important aspect about orthogonality is, operators with 
different transformation properties are necessarily orthogonal. For this reason, the i = 2 
operators are orthogonal to the % = 0, 1 operators, as they transform as (22) and (00) in 
SO (5) respectively. Similarly, e% is orthogonal to ei and ~e~2 as they have different Z7 (10) 
irreps. But ei and as well as f i and g { , have the same transformation properties. 
Can we find a higher symmetry group to distinguish them? We also like to have a group 
such that all the states in d 2 d! belong to a single irrep. 

4-2. The group {7(20), its generators and its subgroup 

The 100 generators of the group {7(10) in terms of the creation and annihilation 
operators are S d v + dftd' v . They can be re-expressed as (Sd)^ SL ' + (d 1 ^ d')^ SL \ where 
S = 0, 1 and L = 0, • • • , 4. This group transforms the states in the configuration d N d' N ' 
among themselves. If we drop the scalar operators for which S = L = 0, we are left 
with the SU (10) generators. The other subgroups that we use, SOs(3) and {7(5), are 
generated by (d)d)^ + (d'td') (10) and (d^d)^ + (d , td') (0L) respectively. 

One can however consider a more general transformation group. If we take the 20 
one-electron states ($ or d'^ acting on the vacuum |0)) as the basis vector space, the 
most general transformation among them constitute the unitary group {7(20), which has 
400 generators of the form d^d v , d'^d' v , d'^d u and d^d' v . These operators preserve electron 
number, but they connect states in different configurations, like d 2 d' with d 3 and so 
forth. The nine operators (e$, <5j and ei) now belong to a single irrep [11 — TT\2o of this 
higher group plus the scalar [0]2o- 

A more elegant formulation is given below. Following Simonis et al (1984) and 
Kaniauskas et al (1987), we introduce the operators d^J^ where t — \ is the isospin 
rank, with d\ sV) = d£ j> , d {t f = d% , d f s = d«»0 an d d\* s = dty . In 

this notation, the 400 generators of {7(20) can be written as w^ TSL ^ = (Sd)^ TSL \ where 
T,S = 0,1 and L = 0, ...,4. One can immediate see that the operators (d^d)( TSL ^ 
with isospin T = are simply the {7(10) operators described above; and those for which 
T = S = correspond to the {7(5) generators. The three components of — ^/l0(d^ d)^ 100 ^ 
are recognized as the pure isospin generators T + ,T and T_. In terms of creation and 
annihilation operators, 

T + = 4^ 

T- = djd„ 

T,= l -(d\d,-d'X) . (15) 

Obviously, the group {7(20) contains SOt(3) x {7(10) as a formal subgroup. 
The subgroup structure {7(20) — > SOt(3) x {7(10) goes exactly the same way as 
{7(10) — ► 505(3) x {7(5). Augmenting equation (||), we can now use the full chain 
of subgroups 

{7(20) D SO T (3) x £7(10) D SO T (3) x 50 5 (3) x {7(5) 
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D S0 T (3) x S0 S (3) x S0(5) D S0 T (3) x S0 5 (3) x S0 L (3) (16) 
^.5. Isospin ranks for states and operators 

Since fermions obey the Pauli exclusion principle, no two electrons occupy the same 
quantum state (spin and orbital quantum numbers m s ,mi, and isospin ^-component 
m t , which is the same as principal quantum number n in a certain sense). Therefore, 
any state with r electrons (d or d') belongs to [l r ] 2 o of Z7(20). Using the branching rule 
given m table |, [l r ] 20 -> r+1 [l r ]io + r_1 [2lC r - 2 )] 10 + • • -, we can quickly write down the 
isospin rank for all the states. Using the same notation as for regular spin, we put the 
multiplicity [T] = 2T + 1 as a superscript preceding the £7(10) label. They are given in 
table [|. Notice that the isospin ranks are uniquely determined by the U(10) labels for 
atomic states. (This is not true for operators, however.) As for the isospin z-component 
Mt, it is simply given by the configuration type. From equation (|T5D , we can see that 
for the configuration d N d' N \ M T = \(N - N'). 

Let us now look at the isospin rank of the Coulomb operators. The operators we 
have considered so far, which do not include CI, all have Mt = 0. This can be verified 
by evaluating the commutator [T ,e ] = for example (using equations (|15D and (|Sc})). 
From now on, we will omit the commutator brackets, knowing that the adjoint action 
(group generators acting on operators) always means taking the commutator; and the 
above equation reads T eo = 0. Since we know that all Coulomb operators belong 
to [11 -TT] 20 and [0] 20 of U(20), the branching rules U(20) -> SO T (3) x 27(10) (from 
table |3|) will give us the possible ranks for the operators. The branching rules uniquely 
determine the isospin rank of the operator e 2 , which belongs to 1 [2 — 2]io- In fact, we 
can easily show that eo and e~\ are also isospin scalars as well. From equation (|TT|), 
we have gj ~ HC^C^d^ d'^d'^d v — d^ d^d^d' u ). Together with equation (|15D , we find 
T + ei = for all i. So, they are all isospin scalars. 

For the other six operators e« and e^, the possible isospin ranks as given by the 
branching rules are 1,5 [0]io and 1,3 ' 5 [11 — 11] 10 . The immediate question that comes to 
mind is whether we can separate the various spin rank components from these operators. 
Repeatedly applying the raising operator T + to them should eliminate the lower rank 
components. One of the results that we find is T_T^e, is proportional to T + ej. This 
tells us that the operators do not have a rank 1 component. So, they belong only to 
1,5 [0]i and ^[ll - TT] 10 . The results are summarized in table 0. 

4-4- Operators with pure isospin rank 

We will now continue our effort to separate the ej (and e«) operators into various isospin 
components. The rank 2 component of is given by ijT^T^ej; as we already know 
there is no spin 1 component, the remaining piece must be the isospin scalar part. From 
equation (fTT|), 

= r (7< + 9i) = a ^{- v dd ■ v d% ■ + ■ v dd> ■ v d/d ■} 

~ k 
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~ E ^c,e(«4^ + 4<H<) • (17) 

Let us take the simplest one, eo = ), as an example. After taking 

four consecutive commutations on eo, we find 

T 2 _Tle = 16e - 4(dJ,4^^ + <!<K<J . (18) 
The first term within the parenthesis is recognized as 2eo (see equation (|3qD). The 
second one is the analogous term for the (i'-electron, which does not contribute in our 
problem as we have only one (f -electron. We would like to absorb the second term into 
our definition of e , and ej in general, and so we write 

Er ( fc ) ( k ) ( fc ) ( fc ) i 

W v dd > ■ : + : vy d , ■ vy d , :} 

~ E ^,(7^(44^^ + «4<) . (19) 

Comparing this result to equation ([FTP , we see that e^ and are almost the same; after 
all, they have the same group transformation properties. Of course, in the complete 
analysis of configurations with at least two d' and two (/-electrons, one should include 
three new operators of the kind d^d^dd — d'^d'^d'd' to account for the interactions between 
the two d'-electrons. 



With our modified definition of eo, the right hand side of equation fll8|) becomes 
16e — 8e , and is a pure isospin rank 2 operator. Repeating the exercise, and keeping 
track of the CG coefficients carefully, one finds that this is true for any i. The 
combination 2ej— e^ has isospin rank 2. We can now separate ej as {|(2ej— ei)+|(e»+ei)}, 
with the first term purely isospin rank 2, and the second one an isospin scalar. One can 
also verify X+(ej + e*) = easily. 

In principle, we can re-tabulate the matrix elements using these six new operators 
Hi + Ci and 2e« — e^ that have simpler transformation properties in place of e« and 
ej. However, we also feel that the relatively simple expressions of e« and e, (as in 
equations (|17D and ( fl9|) ) have their merits. Moreover, the e^ matrix elements are easily 
obtained from the d 2 matrix elements. We would rather leave them untouched. The 
new combinations, with pure isospin ranks, can nevertheless give us a lot of useful 
information and insight into the problem we are facing. 

Recall that we have found some simple proportionality relations in equations ( p.4a| )- 
( |14d) . The underlying reasons are as follows. The combination 2ej — e^, with isospin rank 
2, must have vanishing matrix elements if sandwiched between a pair of 2 [21]io states, 
since two T — ^ states cannot be stretched to T = 2. Therefore, (ej) = |(e.;) for these 
matrix elements, and thus we arrive at relation (|14c ). And the other combination ej + ei, 
being an isospin scalar, must be diagonal in isospin space. So, it vanishes between a pair 
of 4 [111] io and 2 [21]i states. In other words, (e^) = — (ej), as in equation (|14fc|) . From 
our new perspective, the two Wigner-Eckart relations we had before are now interpreted 
as selection rules on the isospin group. 

To explain relation ( |14o| ) requires a more elaborate analysis. First notice that the 
4 [111] io states also appear in the configuration d 3 . In fact, they are almost the same 
state, with the only difference in isospin ^-component. Recall that states in d 3 have 
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isospin Mt = §, while the d?d' states have — \- It is well-known that matrix 
elements can be factored as 

(aTM T \U^\a'T'M' T ) = (-1) T " Mt (_^ t \ ^) (aT\\U^\\a f T') (20) 

where the last factor is the reduced matrix element. So, for the isospin scalar operator 
Hi + ei, we have 

(d 2 d"[lll} l0 p\e t + e t \d 2 d' 4 [lll} l0l ) _ ( | OfW | Of 



-i 



(#4[Hl] lo/5 |e. + e .| d 3 4 [lll]io7) - {-It)!) \-l0 

while for the isospin rank 2 operator 2ej — e», we have 



(dV 4 [lll] 1 o/3|2e t -e i MV 4 [lll] 10 7) _/ § 2 § W | 2§' 
(# 4[ m]lo/ 3| 2 e. _ e .| # 4 [nl]io7) " l_i l) \-lQl 



(21) 



(22) 



2 2/ \ 22. 

But (d 3 4 [lll] 10 /3|ej|d 3 4 [111] 10 7) = in d 3 since there is no d' electron. Therefore 
from these two equations, we can solve for (d 3 4 [lll]i /3|ej|d 3 4 [111] 107) and 
(d 2 d' 4 [lll] 10 /5|e i |rf 2 c/' 4 [lll]io7) in terms of (d 2 d' 4 [lll] 10 /3|e l |rf 2 d / 4 [lll]io7) which are 
assumed known; and the answer is 

(d 2 d' [lll] 10 /3|ei|dV [lll]ioT> = 2 (rfV [Hl]io/9|e^V [lll]i 7> (23a) 
(d 3 [lll] 10 /?|e 4 |c/ 3 [lll]io7) = 3(dV [lll]io/3|ei|dV [111] 107) • (236) 

The first one is of course relation (|14a|) that we are aiming at; the second one is a bonus, 
that relates the d 2 d' matrix elements to the ones in d 3 . Of course, the d 3 matrix elements 
are also available in the literature (see Nielson and Koster 1963 for instance), so the 
second equation can be served as a consistency check. 

4-5. Configuration mixing via Coulomb interaction 

The isospin arguments used above can relate the d 3 Coulomb matrix elements to those in 
d 2 d'. In the same way, we can obtain the mixed matrix elements between configurations 
d 3 and d 2 d' . The first point to note is that the Coulomb operators which connect d 3 with 
d 2 d' necessarily have Mt = 1; so no isospin scalar Coulomb operator can be responsible 
for CI. Let us write (Ai)^ 2 J as the isospin rank 2 Coulomb operator, such that 

G£) { ? = y|(2e, - e,) 

/2 11 
3 C^C^S^d'^dy + d\d^d! v - ^d\d\d^d v - -dJdJ4<,) 

= C 1 pyCri^d^d^dpdv — d^d^d^d^ 

(A) ( 2 = -J2 C ^ C nd d l d l d W ■ 
With the notation developed before, the above operators are simply the components of 
the tensor operator 10 (w/ 10 ^ • u!( 10fc ))( 200 ). Parallel to equations (21) and (|22]), we have 



{d* *[Ul] w l3\(AfM* Wiol) ,_ f I 2|W § 2 j\ _ 1 



(d 2 d' 4 [iii] 10 /3i(^) (2 oi^ /4 [iii]io7) v-l 1 !/ v-l°|y V2 



(24a) 
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Table 7. Eigenvalues of the U(10) scalar operators. 

operator eigenvalue 

e |iV(iV - 1) + ±N'{N' - 1) 

7o NN' 

g : T+T_ : = T 2 - T§ + T - N 



= T 2 — + l) + TViV' 



i {T 2 _ ( iv±w: ) (iv±A L : + 1)} + ^ 



2 ' V 2 

|{_ T a + (%Ci)(%Cl+ i, J 
±{T 2 + 3( 
2e - e T 2 - 3T 2 



eo + eo I{T 2 + 3(^)(™i-l)} 



<d 3 4 [lll] 1 o J g|(^) ( ?|d 2 d / 2 [21] 10 ^> = _/ | 2|W § 2± 

(d 2 d' 4 [lll]io/3|(^) ( SVd ,2 [21]io7> " V-l 1 l7 V-|°l 

The other possible configuration mixing Coulomb operators (My = 1) are given by 
( w (iofc) . w (iofc))(ioo) and ^)(i) _ ( w (oofc) . w (iofc))(ioo)_ The firgt one ig in fact identically 

zero when normal ordered. The second one with My — are the isospin rank one 
operators (i = 0, 1,2) that we briefly mentioned in section [4.4| . They are of the form 
d^d^dd — d'^d^d'd', as opposed to the operators in ([19]). In the configurations with 
no more than one d' electron, they have the same matrix elements as the e^. So in this 
manner, we can determine all the configuration-interaction Coulomb matrix elements 
between d 3 and d 2 d' from the known ones in d 2 d'. We can similarly write down the other 
matrix elements between d 3 and dd' 2 , d 2 d' and dd' 2 and so forth. 

Before we close this section, we should mention the relevance of Brillouin's theorem 
(see Bauche and Klapisch 1972, and Godefroid et al 1987) in our analysis. The theorem 
says, the Hartree-Fock (HF) solution ^ hf of the configuration l N has vanishing matrix 
elements with a class of states in Z^" 1 /'. They are the states in which the /' electron 
is coupled to the l N_1 state via the ordinary fractional parentage coefficients as in l N . 
That is to say, states in l N ~ x V and I with the same angular form will not mix via 
Coulomb interaction, if the HF solution for the l N configuration is used. In our case, 
(d 3 4 [lll] 10 /9|Coulomb|(i 2 (i / 4 [lll]i /9) will have a vanishing radial integral. That does 
not make our work less useful, however. In order to apply Brillouin's theorem in our 
situation, one has to first find out the HF solution ^hf for each single LS term in 
the configuration I , which might be straightforward, but definitely not easy. Only 
then, equation ( p4o|) will rendered irrelevant, as Brillouin's theorem predicts zero radial 
integrals in those cases. All other results are otherwise unaffected. 

4-6. Eigenvalues for the scalar operators 

From table ^, there are a few complete scalar operators. We can now find closed- 
form expressions for each £7(10) scalar operator; that is, the i = operators. Let us 
start from the simplest operator eo- From equations (|8~aD and (|19D , we can see that 
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the eigenvalue for eo is simply given by ^N(N — 1) + ^N'(N' — 1). For f , from 
equations (ffij), it is equals to NN'. For the next one, g , from equations ( |8c|) and (p~5|) , 
we find g = : T + T_ : = T + T_ — AT. None of these three operators e , f and g is an 
isospin scalar. The other two operators eo and e~o are just linear combinations of f and 
g . The results are summarized in the table [?]. Note that the operator e~o is an isospin 
scalar, so is the combination cq + eo- Their expressions from table [7] certainly verify that 
fact; as T 2 is a scalar in SOt{3) x Z7(10) and the total electron number = N + N' 
is, in fact, a scalar in £7(20) and hence scalar in all its subgroups. 



5. The operator sum = + f i = + + 

Recall that at the end of section |3.2| , we mentioned a simple result on the operators 
gj = ej + f { . In terms of the new operators, the sum is e; = + + e^; whose matrix 
elements are included in table [5[ As a brief summary of our work up to this point, the e^ 
matrix elements of d 2 d' are obtained easily from the ones in the configuration d 2 . The e, 
matrix elements can be related to e, by just three constants using the isospin structure 
as in (|14a| ) - ( |14dj . For the last operator e iy selection rules on 50^(3) x Z7(10) give us a 



lot of vanishing matrix elements (equations ( |1 3 oj) and ( |135|) ); the rest involving a pair of 



[21] 10 states are yet to be found. A thorough understanding on the simple result for the 
operator e, will definitely help to accomplish our plan for finding the remaining matrix 
elements. 

Let us look at the three different cases separately. For i = 0, e always equals 3 in 
d 2 d! . From table 0, one can see that the operator sum equals ^Nt(Nt — 1), which is 
scalar in U(20). Since all the states belong to the same U(20) multiplet [lll] 2 o, they 
must have the same matrix element. 

The next one e\ is also diagonal; with matrix elements 0,3,6,7 and 13 only, which are 
determined by the U(5) and 5*0(5) irreps. To explain this, we write down the operator 
in terms of = v^J + v^\, = 2w^ 00k \ Using the extended version of equations (ff^ ) 
(recall that we have added the identical terms involving the cf-electron to the operators 
Ci in section |4^4j), the operator e\ takes the simple form 

Z . „(o) . „(o) . + :v m. v m. + :„(*).„(*); 

2 

which is reminiscent of the (quadratic) Casimir operators for the groups SU(2l + 1) and 
50(2/ + 1); they are given by C(SU 2l+1 ) = • «(*) and C(S0 2l+1 ) = £ «(*> ■ 

k>0 odd k 

respectively. For the irrep [Ai, A2, ■ • •] 2Z+1 of U{21 + 1), the eigenvalue of C(SU2i+\) is 
(see Judd 1998, §5.9) 

21+1 2 

X)A i (A i + 2 + 2Z-2i)- T -r (25) 

i=l rJ 
where n — X) K- One can check that [Ai, A 2 , • • •] 2i+i an d [Ai + a, \ 2 + a, ■ ■ -]2z+i both 
have the same eigenvalue as they should since they possess the same 577(2/ + 1) content. 



Coulomb and spin- orbit interaction matrix elements in d 2 a" configuration 



21 



The eigenvalue of C(S , 2 ;+i) on the irrep (wi, w 2 , • • •) of £0(2/ + 1) is 
1 1 

-X>i(wi + l + 2Z-2i) . (26) 
1 i=i 

We can remove the normal ordering using the relation : ■ : = ■ «W — M Nt, 
where I = 2 in our case. Putting these together with ■ = N^/[l], we find 

ei = ^ : iV| : + : C(£{7 5 ) - C(£0 5 ) : 

= ^iV| - \n t + C(£[/ 5 ) - C(50 5 ) . 

This explains why e.\ is diagonal, and gives the correct eigenvalues. As we know (ei) and 
(ei) already, this simple result can also be used to determine (ei) = in equation (|13a|) , 
where vanishing matrix elements cannot be explained by Z/(10) selection rule. 
Finally, the operator e 2 is almost diagonal, with the three exceptions 

(d 2 d' ^U[21}(21)^D\e 2 \d 2 d' ^U[21](10)^D) = 6^21 . 

Since the operators e 2 + e 2 and e 2 are both isospin scalars, clearly the sum must also 
be diagonal in isospin space. With a little hindsight, we can further conclude that it is 
diagonal in the {7(5) space as well; so that the [21] 5 states will not mix with [111] 5 or 
[3] 5 states^. The above three exceptions are the only possible off-diagonal entries we can 
have. Furthermore, the matrix elements for a set of L states do not depend on the spin, 
isospin or U(10) irreps; they depend only on their respective £0(5) irreps. All these 
can be explained using the spin-isospin supermultiplet group £{7(4) of Wigner (1937). 

To explain these, let us first introduce an alternative branching scheme {7(20) D 
£O s (3) x Z7'(10) D £O s (3) x £0 T (3) x {7(5) as oppose to the one in flU). The U'(10) 
group acts on the isospin-orbital space, which is analogous but different from the spin- 
orbital {7(10) group that we had before; however, the {7(5) group is the same in both 
schemes. To better display the symmetry, in the original scheme, we label a state with 
its spin superscript [£] put before the {7(5) label. For example, we now write the first 
of the doublet F states and the last of the quartet F states as | 4 [111] 10 2 [21] 5 (21)F) 
and | 2 [21] 10 4 [21] 5 (21)F); which become | 2 [21]' 10 4 [21] 5 (21)F) and | 4 [111]' 10 2 [21] 5 (21)F) 
respectively in the new scheme. With a mere spin-isospin exchange, the doublet F 
state switches with the quartet F state; on the other hand, the operator e 2 is invariant. 
This explain why the two states have the same matrix element. But there is yet another 
doublet F state, | 2 [21] 10 2 [21] 5 (21)F), that has the same diagonal matrix element. To 
accommodate this, we should promote the idea of spin-isospin exchange to a more 
general transformation in the spin-isospin space. This leads us to introduce the £{7(4) 
supermultiplet group due to Wigner (1937). The three F states now fall into a single 
supermultiplet [21] 4 ® [21] 5 of £{7(4) x {7(5) C {7(20)|. The other possible multiplets 
in d 2 d' are [111]4 <8> [3] 5 and [3] 4 ® [llljs. The operator e 2 , being a £{7(4) scalar, is 

* The operator is not a {7(5) scalar however; U(b) scalars do not give (22) of SO(5). 

t) We use the simpler U(4) irreps rather than SU(4) irreps throughout the paper, which should not 

affect the validity of our arguments and results. 
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diagonal in the SU (4) space; which appears as if it is diagonal in the orbital £7(5) space 
since each £7(5) irrep is paired with a unique SU (4) irrep in the problem on hand. This 
also explain why the matrix elements are independent of spin or isospin ranks, as they 
belong to the same multiplet in SU(4) x SO (5). 

As explained before, (d 2 d' 2 [21]i p\e2\d 2 d /2 [21] W 'y) are the only matrix elements for 
which we have not yet found a simple way to deduce the values. With this latest result, 
we can obtain those matrix elements easily if /3, 7 do not contain the irreps [3]s(30) of 
£7 (5) and SO (5). To fill in the last piece of puzzle, we make use of the following operator 
introduced by Judd (1998, p 222): 

n' = 7v^ • «W - 3v^ . „(3) = -3C(S0 5 ) + L 2 . 

It transforms as (22) in 5*0(5), just like the i — 2 operators; and its eigenvalue is easy 
to compute. So for the | 2 [21]i [3] 5 (30) 2 SG.H7) states, since (30) x (30) contains (22) 
only once, we can apply the Wigner-Eckart theorem and find 

( 2 [21][3](30) 2 L|e 2 | 2 [21][3](30) 2 L> = ( 2 [21][3](30) 2 L|fif[21][3](30) 2 L) . 

In other words, we only need to calculate the matrix element for one of the four terms 
SGHI to determine the proportionality constant (which turns out to be unity). The 
other three can be determined from the above relation easily. Usually, the calculation 
on the fully stretched state ( 2 J) is reasonably easy. In the end, we find that ( 2 /|/ 2 | 2 -0 
is the only matrix element we need to calculate. All other d 2 d' matrix elements are 
related to e, ones via some simple arithmetic relations. This is a truly surprising result 
from the present analysis. 



6. Spin-Orbit interaction 

From a group-theoretical point of view, the spin-orbit interaction is much simpler than 
the Coulomb interaction. The interaction Hamiltonian, in second quantized form, is: 

tt /- (11)0 , J- (11)0 

H 2 = (m w dd + Cd'd' w d ,J 

where the (dd and (d'd. 1 are the corresponding radial integrals; the symbol k following 
the ranks (SL) for the coupled tensor w dd K = — (dftd)( SL ^ K is the rank to which S and 
L coupled (see Judd 1967). The radial integral (dd is related to the classic parameter ( 
of Condon and Shortley (1953) by (dd = —VTE(- In the same manner, we can include 
CI into our analysis by introducing the perturbing Hamiltonian 

rrmix _ f „ (11)0 > (11)0 

n 2 - Ud' w dd , + Qd'd w d , d . 

The operators and w d ) d ^° are simply hermitian conjugate to each other. Each 

of the above four operators belongs to [1 — T^cb the adjoint representation of U(20); 
since the nine components of are the SU (20) group generators (see section |4.2[ ), 
the coupled scalar operator must also belong to the adjoint representation. Using the 
branching rule £7(20) -> SO T (3) x £7(10) given in table |, [1 - Y] 20 gives x ' 3 [l - T] w + 
1,3 [0]i . Only the irreps 1,3 [1 — l]i can eventually give us a S = L = 1 ( 3 P) tensor in 
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the ordinary spin and orbital spaces, to which the spin-orbit operators belong; hence the 
Z7" (10) scalar irreps can be dropped. Then from C/ (10) — > 50,g(3) x U(5), [1 — l]io gives 
us 1,3 [1 — 1)5 + 1,3 [0]s; once again, the U(5) scalar irreps can be dropped. And this time 
we only want the spin 1 part, so 3 [1 — Tjs is the only permissible label for the spin-orbit 
operators on the 5*05(3) x U(5) level. 

Finally, on the 50(5) level, using the branching rules £7(5) — > 50(5), [1 — 1] 5 gives 
us (11) + (20) of 50(5), but (20) contains only 5, D and G terms in 50 L (3) while we 
are looking for a P term. So, we can conclude that the spin orbit operators belong to 
[1 — Tj2o 1,3 [l — T] 10 [1 — T]5(11) 3 P in our branching scheme. 

6.1. Relations between various spin-orbit matrix elements on the C/ (10) level 

Using the Wigner-Eckart theorem, we can take advantage of the known spin-orbit 
interaction matrix elements of the configuration d 3 (see Nielson and Koster 1963, 
for instance), in which all the states belong to 4 [111] 10 of 50^(3) x [7(10). Since 
[111] 10 x [lll]io contains the adjoint irrep [1 — l]i once, on the £7(10) level we can 
write down several direct proportionalities. With the numerical constants, they run 

(d 3 IlllMII^N 3 4 [Hl]io7) = 3(dV ^lll] 10 (3\\w^\\d 2 d' 4 [111] 10 7) (27a) 

= 4 [lll] 10 /3||^ 1) ||rfV 4 [111] 10 7) (276) 

= ^{d 3 4 [lll] 1 o/3|ki 1 i ) ||rfV 4 [111] 10 7) • (27c) 

Since these involve reduced matrix elements, no components of u/ 11 ) or magnetic 
quantum numbers Ms and Ml appear; hence the extra information k — is rendered 
irrelevant. So, with just three constants, we can determine all the spin-orbit matrix 
elements between a pair of 4 [111] 10 states in d 2 d', as well as the CI matrix elements 
between d 2 d' and d 3 . 

For matrix elements involving the 2 [21]i states, there are no such corresponding 
states in the d 3 configuration. But still, for matrix elements between the 4 [lll]i and 
2 [21] 10 states, since [llljio x [21]i contains the adjoint only once, we can relate the 
reduced matrix elements by direct proportionalities, which turn out to be 

(d 2 d' '[lll] w (3\\w^\\d 2 d' 2 [21] lo7 ) = ~(d 2 d' 4 [lll] 10 /?||^i ) ||dV 2 [21] lo7 ) (27d) 

= ~^(d 3 *[lll] w p\\w^\\d 2 d' 2 [21] lo7 ) (27e) 

Lastly, for matrix elements between a pair of 2 [21] 10 states, the simple version of the 
Wigner-Eckart theorem does not apply, since [2T] 10 x [21] 10 contains the adjoint [1 — T]io 
twice. The above equations are all we can get from the Wigner-Eckart theorem on the 
£7(10) level. One direction to proceed is to seek for similar relations on the U(5) and 
50(5) levels, so that most of the remaining matrix elements can eventually be related 
to the ones in d 3 . We can also look at the use of isospin in the spin-orbit analysis. 
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6.2. The use of is o spin in spin-orbit interaction 

In order to take advantage of the underlying isospin structure in the spin-orbit 
interaction, we should first decompose the spin-orbit interaction operators into 
components with definite isospin ranks. The four operators w]^\ , w^' and 
belong to 1,3 [1 — T\ w in SOr(3) x U(10). Among them, v)^, and w^j has M T = 1 and 
— 1 respectively; so they must belong solely to 3 [1 — T]io- For the other two operators, 
and WjijJ, with Mt = 0, one can check that + w^J is an isospin scalar 

(i.e. commutes with T+/_); and w^P — w d }^} has an isospin rank 1. More precisely, 
—w dd ]\ ^( w dd ~ w &d') anc ^ w d/d are ^ ne three components M T = 0, ±1 of the operator 
3<u;( m ) = 3(d) dY 111 ^ ; while —y/^w^- 011 ^ = (w^ + w^'d') i s ^ ne isospin scalar operator. 

We can now easily deduce equations (|27a|) to (|27^) from isospin arguments. First 
of all, the scalar operator just mentioned is diagonal in isospin space, therefore we have 

{a?d ^lllj^ll^ + w d ^,\\d 2 d 2 [21] lo7 ) = 



which proves relation (]27d|) . Then, using equation ( pL|) to remove the dependence 
of the two matrix elements 

(rf 3 %ll\ w (3\\w^)\\d?d> 2 [21] lo7 ) {= --^(Illll^lllS^^IH 2 ^]^)^ 
(d 2 d' "[ni] 10 ^(w^ -w^})\\d 2 d' 2 [21] lo7 ) 



we can easily get relation ([27 Similarly, we can express 
(d 3 '[in] 10 ^\\±(w^+w^)\\d 3 4 [111] 10 7) 

(d 2 d' Illllio^H-^^ +w^)\\d 2 d' 4 [111] 10 7) 

(d 3 *[in] 10 f3\\-L(wW -wW)\\d* 4 [lll] lo7 ) 

(d'd' Illllio^H-^^ -w^})\\d 2 d' 4 [111] 10 7) 

(d 34 [lll] 1 o/5||^ ) ||d 2 d' 4 [lll] 10 7) 

in terms of their matrix elements reduced with respect to S, L and T. The first four 
matrix elements, together with the fact that {w d ) d ,) = in d 3 , give equations ( |27q| ) 
and (|276| ). The last one gives the remaining relation 



7. Concluding remarks 

We started with this relatively simple configuration d 2 d' as the subject of the present 
work. The ideas and techniques introduced in this paper can certainly be applied to 
other configurations like d N d' and d N d' N ' . Of course, it is a different question whether 
these configurations are experimentally observed, or simply of theoretical interest. 
Another obvious extension is to apply it to the /-shell electrons. The /-shell is known 
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to have a richer group structure. We are aware that similar analysis on the /-shell might 
reveal more surprises calling for explanations. 

Another direction for further research is to introduce the idea of quasispin into 
the analysis. Quasispin generators are of the type (cftcfty 00 ^ and (dd)^ in the case 
of ordinary ci-electrons. These operators, together with (<fl~<i)( 00 ) + (dtfi)(°°\ generate 
the SOq(3) quasispin group. But in the case where the electrons have an extra 
spin structure, namely the isospin, it is known that the analogous quasispin group 
is qualitatively different (see Flowers and Szpikowski 1964, Feng and Judd 1982). The 
full quasispin group is 5*0(8) rather than merely SO (3). This extra group structure 
may gives us some new selection rules and proportionality relations, as we know it does 
in the case of configurations of equivalent electrons. 
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